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Abstract 

It has been widely known that bino-like dark matter in the supersymmetric (SUSY) the¬ 
ories in general suffers from over-production. The situation can be drastically improved if 
gluinos have a mass slightly heavier than the bino dark matter as they reduce the dark matter 
abundance through coannihilation. In this work, we consider such a bino-gluino coannihilation 
scenario in high-scale SUSY models, which can be actually realized when the squark-mass scale 
is less than 100-1000 TeV. We study the prospects for exploring this bino-gluino coannihilation 
scenario at the LHC. We show that the searches for long-lived colored particles with displaced 
vertices or large energy loss offer a strong tool to test this scenario in collider experiments. 



1 Introduction 


The first stage of the LHC running has pointed a possible direction for the actual realization of 
the supersymmetric (SUSY) Standard Model (SM). First and foremost, the observed SM-like Higgs 
boson [1] with a mass of about 125 GeV [2] implies that the mass scale of SUSY particles is higher 
than the electroweak scale; the radiative corrections by stops easily lift up the Higgs mass from 
the tree-level value predicted to be less than the Z-boson mass in the minimal SUSY SM [3], if 
the stop masses are far above the electroweak scale [4, 5]. This is in fact consistent with lack 
of any evidence in the SUSY searches so far [6, 7]. A relatively high SUSY breaking scale offers 
further advantages for SUSY SMs. For instance, heavy masses of SUSY particles suppress the flavor 
changing neutral current processes as well as the electric dipole moments of the SM particles [8, 9], 
which are stringently constrained by the low-energy precision experiments. Moreover, such heavy 
SUSY particles reduce the proton decay rate via the color-triplet Higgs exchange [10] and make 
the simplest version of the SUSY grand unihcation model [11] viable. In cosmology, the gravitino 
problem is evaded when the gravitino mass is high enough [12]. These attractive points stimulate 
quite a few studies of high-scale SUSY models [13-21]. 

An order parameter of SUSY breaking is the gravitino mass m 3 / 2 . If the SUSY breaking effects are 
transmitted to the visible sector via the gravitational interactions (or other interactions suppressed 
by some high-scale cutoff such as the Planck scale), then the soft SUSY-breaking scalar masses 
are induced with their size being 0 ( 7713 / 2 ). In this case, the scalar SUSY particles typically have 
masses of the order of 7773 / 2 ; from now on, we express the typical masses of these scalar particles 
by m ~ 7773 / 2 . The masses of the fermionic SUSY particles (gauginos and Higgsinos) are, on the 
other hand, dependent on models, since their mass terms can be suppressed if there exist additional 
symmetries. For example, the gaugino masses become much smaller than the gravitino mass if the 
SUSY breaking helds are charged under some symmetry. In this case, these masses are generated 
by quantum effects, such as anomaly mediation contribution [22, 23] and threshold corrections at 
the SUSY breaking scale [22, 24]. They are also affected by the presence of extra particles [25]. 
Moreover, the Higgsino mass can be suppressed by, e.g., the Peccei-Quinn symmetry [26] and be 
much lighter than 7773/2 and m. See for instance Refs. [27, 28] for a concrete realization of light 
Higgsinos. 

Possible deviation of the masses of the fermionic SUSY partners from 7773/2 and rh gives additional 
benehts to SUSY SMs. Firstly, if gauginos he around 0{1) TeV, gauge coupling unihcation is realized 
with great precision [29] even when the scalar mass scale m is much higher than the electroweak 
scale. Secondly, the neutral components of these fermions, the neutral bino, wino, and Higgsino, 
can be a candidate for dark matter (DM) in the Universe. Among them, the neutral wino is one of 
the most promising candidates since the anomaly mediation mechanism naturally makes the wino 
be the lightest SUSY particle (LSP). Its thermal relic abundance actually explains the observed 
DM density if the wino mass is around 3 TeV [30]. Currently the mass of the wino LSP is 
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restricted by the direct search at the LHC as > 270 GeV [31]. The wino DM scenario is also 
being constrained by the indirect DM searches nsing gamma rays [32, 33]. These experiments, as 
well as the DM direct detection experiments [34], can probe this scenario in fntnre. Higgsino DM 
with a mass of ~ 1 TeV can also acconnt for the observed DM density [35] . For the recent stndy of 
the phenomenology and fntnre prospects for this Higgsino DM scenario, see Ref. [36] and references 
therein. 

The last possibility is bino DM. If the scalar SUSY particles and Higgsino are signihcantly 
heavy, bino DM is usnally over-prodnced as the interactions of bino with the SM sector tend to be 
snppressed. To avoid the over-prodnction and get correct dark matter abnndance, we need some 
exceptional mechanism to rednce the bino abnndance, such as coannihilation and Higgs funnel [37]. 
If the Higgsino mass is heavier than 0(10) TeV, the remaining possibility is the coannihilation. 
In this case, its thermal relic agrees to the observed value if there exist some particles degenerate 
with the bino DM in mass. In fact, as shown in Refs. [38-43], bino DM can explain the correct 
DM density if wino or gluino has a mass slightly above the bino mass. After all, there are various 
options for DM candidates in the high-scale SUSY scenario, and therefore it is quite important to 
experimentally examine each possibility. 

Among the possibilities mentioned above, the collider testability of the bino-gluino coannihilation 
is expected to be the most promising since this case requires light gluinos. As we shall see below, we 
expect an 0{1) TeV gluino mass in this case, which can be within the reach of the LHC. This could 
be compared to other DM scenarios in high-scale SUSY models; for instance, if the gaugino masses 
follow the spectrum predicted by the anomaly mediation, wino is the LSP and it becomes the main 
component of DM if it has a mass of 3 TeV, as mentioned above. In this case, the gluino mass 
is predicted to be 0(10) TeV, which is of course far above the possible reach of the LHC. In this 
sense, it could be much easier to look for gluinos in the bino-gluino scenario than other cases. This 
naive expectation, however, turns out to be questionable. The bino-gluino coannihilation scenario 
requires that the mass difference between bino and gluino, AM, be AM < 100 GeV. Such small 
mass difference results in soft jet emissions, which make it extremely challenging to detect the signal 
of gluino production. For this reason, previous studies have concluded that it is difficult to probe this 
bino-gluino coannihilation scenario at the LHC if the DM mass is heavier than 1 TeV [42, 44, 45]. 

In this work, we show that this small mass difference actually helps us to probe the bino-gluino 
coannihilation. When AM < 100 GeV and the sfermion masses are much heavier than the gaugino 
masses, the lifetime of gluinos Tg can be long enough to distinguish its decay signal from that of 
prompt decay. As will be shown below, we expect its decay length to be CTg > 0{1) mm when the 
sfermion masses are 0(100) TeV. A decay length of this order is in fact the main target of searches 
for long-lived colored particles with displaced vertices (DVs) [46] and large energy loss [47]. We 
will hnd that this search technique indeed gives a stringent limit on the bino-gluino coannihilation 
region, and probe wide range of the parameter space in future experiments. 

This paper is organized as follows. In the next section, we consider the bino-gluino coannihilation 
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scenario and show the parameter region which accomplishes the correct DM density. The lifetime of 
gluino predicted in this parameter region is given in Sec. 3. Then, in Sec. 4, we discuss the strategy 
of the long-lived gluino searches at the LHC, and present the current constraint and future prospects 
for the bino-gluino coannihilation scenario. Finally, Sec. 5 is devoted to conclusion and discussion. 

2 Bino-gluino coannihilation 

To begin with, let us discuss the bino-gluino coannihilation scenario [39-43] to clarify the target 
parameter space we consider in the following analysis. Throughout this paper, bino is assumed to be 
the LSP and be the DM in the Universe. We consider the case where the bino-gluino coannihilation 
is effective so that the thermal relic abundance of the bino LSP is consistent with the observed DM 
density = 0.12. Thus, bino and gluino should be degenerate in mass, i.e., AM = Mg — < 

100 GeV, with Mg and M^ being the gluino and bino masses, respectively. We further assume that 
the typical mass of scalar SUSY particles, m, as well as the Higgsino mass Mjj, is as high as the 
gravitino mass 7713 / 2 . This setup is realized with a generic Kahler potential. The gaugino masses 
are supposed to be suppressed by a loop factor compared with 777 . 3 / 2 , which occurs when the SUSY 
breaking superhelds are non-singlet. Namely, we require M^ ~ Mg m Mfj ~ 7773/2 in what 
follows. Moreover, we assume the wino is heavy enough not to contribute to the coannihilation 
process. It turns out that such a mass spectrum can be in fact realized in the high-scale SUSY 
models [25, 42]. We will see below that the scalar mass scale m gives the signihcant effects on the 
determination of the bino DM abundance.^ 

The relevant annihilation processes to the computation of the thermal relic abundance are the 
self-annihilation and coannihilation of bino and gluinos. Among them, gluino self-annihilation is the 
most effective because of the strong interaction, and this plays the dominant role in the determination 
of the bino relic abundance. The bino self-annihilation and bino-gluino annihilation are much smaller 
than the gluino self-annihilation, since these cross sections are suppressed by heavy Higgsino and 
sfermion masses. Hence, these annihilation processes scarcely affect the following calculation. 

An important caveat here is that the bino-gluino coannihilation does not work efficiently with¬ 
out chemical equilibrium between bino and gluinos [43, 48]. Therefore we should require that the 
transition rate between them should be fast enough compared to the Hubble expansion rate. The 
transition rate is, however, again suppressed by heavy squark masses. Thus, we obtain an upper 
bound on in by imposing the above condition. The transition rate of bino into gluino via quark 
scattering, T{Bq —)■ gq), is estimated by the product of the corresponding scattering cross section, 
(T{Bq —)■ gq), and the number density of initial state quarks, Uq. The former is approximately given 
by (T{Bq —)• 'gq) ~ T‘^/ih^ with T being the temperature of the Universe, while the latter is Uq ~ T3 
since quarks are relativistic when the transition process is active. Consequently, the transition rate 

^While completing this manuscript, we received Ref. [43], which also discusses the squark mass effects in the gluino 
coannihilation scenario. 
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is given by 


2^5 

V{Bq ^gq)^—^ 




( 1 ) 


On the other hand, the Hubble rate H goes like H ~ T^/Mpi with Mp\ the Planck scale in the 
radiation dominated epoch. In order to sufficiently reduce the bino density through coannihilation, 
the condition T{Bq —)■ 'gq) 3> H should be satisfied until the bino DM decouples from thermal bath 
at the freeze-out temperature T/ ~ M^/20. This reads 
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which then gives an upper bound on the scalar mass scale m. Here Xf = M^/Tf ~ 20. Numerically, 
we have 

/ M- \ 4 

^ < 250 X 

We find that when the DM mass is 0{1) TeV the upper bound on the scalar mass scale lies around 
0(10(2-3)) TeV; indeed, many high-scale SUSY models [13-21] predict the SUSY breaking scale to 
be this order, with which the 125 GeV Higgs mass is naturally accounted for. Therefore, it is quite 
important to take into account the constraint on m when we discuss the bino-gluino annihilation in 
the high-scale SUSY scenario. 

To make the above discussion more accurately, we perform the numerical computation by solving 
the Boltzmann equation to obtain the bino-gluino conversion rate and the resultant relic abundance. 
First, in Fig. 1, we show the ratio of the bino-gluino conversion rate Fwith respect to the 
Hubble rate H as functions of Mj^/T. Here, we set M^ = 1.5 TeV, AM = 50 GeV and m = 
(200,300,400) TeV in the red solid lines; Mg = 0.5 TeV, fh = 300 TeV and AM = (50,100,200) GeV 
in the blue dashed lines; = 3 TeV, in = 100 TeV and AM = (50,100, 200) GeV in the green 

dotted lines. All of the squark masses are assumed to be equal to the universal mass in. When 
we evaluate the transition cross sections and (inverse) decay rate of gluino and bino, we use the 
effective theoretical approach to properly deal with sizable quantum corrections resulting from large 
difference between the gluino and squark mass scales; we first integrate out squarks to obtain a set 
of dimension-six operators which involve quarks, bino and gluino, and then evolve these operators 
down to the gluino mass scale by using the renormalization group equations, which results in a 
several tens percent enhancement of the transition rate, compared to the tree level calculation 
[49-51]. The loop-induced dimension-five dipole operator (gluon-bino-gluino) is found to be quite 
suppressed and thus its contribution is negligible in the present analysis. In addition, we include the 
so-called Sommerfeld effects [52] on the gluino annihilation. On top of that, p-wave contribution, 
finite-temperature effects, the scale dependence of the strong coupling constant in the QGD potential 
[41], possible ambiguity in the initial state color arrangement^ due to thermal effects [53], and the 

^In our computation, we assume that the initial state gluinos have a definite color configuration, not thermal 
averaged one. 
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Figure 1: Ratio of the bino-gluino conversion rate to the Hubble rate as functions of M^/T. We set 
Mg = 1.5 TeV, AM = 50 GeV and m = (200,300,400) TeV in the red solid lines; Mg = 0.5 TeV, 
m = 300 TeV and AM = (50,100, 200) GeV in blue and dashed lines and rrig = 3 TeV, m = 100 
TeV and AM = (50,100, 200) GeV in green and dotted lines. 

bound-state effects on a pair of gluinos [43] may change the results by a factor of (9(10)%. The above 
figure shows that the conversion rate decreases as rh or AM is taken to be larger. In particular, if 
the squark mass scale rh is several hundred of TeV with the DM mass being a relatively small, then 
the condition T g^- ^ H does not hold any more when the DM abundance freezes out. 

In Fig. 2, we plot on the Mg — rh plane the mass difference AM with which the thermal relic 
abundance of bino DM explains the observed DM density = 0.12. In the red shaded region, 

the squark mass is too heavy for the coannihilation process to work well and therefore the DM is 
overproduced. We will discuss how to probe the parameter space shown in Fig. 2 at the LHG in the 
subsequent section. 

3 Gluino lifetime 

Next, we study the lifetime of gluino, which plays a crucial role in the discussion of the testability 
of the bino-gluino coannihilation scenario at the LHG in the following section. As mentioned in 
the Introduction, in this scenario, a relatively light gluino mass is expected. Thus, the gluino pair 
production is suitable target for the hadron collider experiments like the LHG in this case. After the 
pair production, a gluino decays into a bino, a quark, and an anti-quark through the squark-exchange 
processes [49-51, 54]. When the gluino is degenerate with the bino in mass, which is required in the 
bino-gluino coannihilation scenario, the decay length of the gluino, CTg, is approximately given as 
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Figure 2: Contour for the mass difference AM which makes the thermal relic abundance of bino 
DM equal to the observed DM density Q-ouh^ = 0.12. In the red shaded region the bino DM is 
overproduced due to failure of bino-gluino coannihilation. 


follows: 


CTg = 0(1) X 


AM 

100 GeV 


-5 


m 


100 TeV 


cm . 


(4) 


From this equation, we see that the decay length gets longer as the mass difference AM is taken 
to be smaller or the scalar mass scale rh is set to be larger. Therefore, we expect a relatively long 
decay length when the bino-gluino coannihilation is achieved in the high-scale SUSY scenario. 

To illustrate the gluino decay length corresponding to the bino-gluino coannihilation region, in 
Fig. 3, we plot contours of the gluino decay length in colored lines with the squark masses set to 
be m = 100 TeV, which we denote by on the M^-AM plane. We also show the mass 

difference AM with which the thermal relic of the bino DM agrees to = 0.12; the black solid 

line shows the case where the bino-gluino chemical equilibrium is assumed, while the other black 
lines represent the cases of fh = 100, 300 and 500 TeV. To avoid overproduction, AM should be 
below these lines. From Fig. 3, we hnd that the gluino decay length is scarcely dependent on the 
bino mass, which has been already shown in Eq. (4) implicitly. We have CTg > 0{1) mm where 
the thermal relic abundance of the bino DM explains the observed DM density. This is a crucial 
observation for the strategy of exploring the bino-gluino coannihilation region at the LHC. 
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Figure 3: Decay length of the gluino with the squark mass m = 100 TeV in colored (almost 

horizontal) lines. Mass difference AM with which the thermal relic of the bino DM agrees to 
= 0.12 is also shown in the black solid line for the case in which the bino-gluino chemical 
equilibrium is assumed, while the cases for m = 100, 300 and 500 TeV are given in the other black 
lines. 


4 LHC search 

If gluino decays promptly and the bino and gluino masses are almost degenerate, it is quite hard to 
search for the gluino at the LHC, since the small mass difference makes the missing energy and jet 
activities tiny. Currently the ATLAS and CMS collaborations have put limits on such a degenerate 
neutralino, he., bino in our case, with a mass of around 600 GeV [6, 7]. The bounds are expected 
to reach ~ 1200 GeV with the integrated luminosity of 300 fb~^ at the 14 TeV LHC [55]. 

These limits are in fact drastically improved once we consider the fact that in the case of the 
bino-gluino coannihilation scenario, the gluino lifetime is as long as CTg > 0{1) mm, as we have seen 
in the previous section. Such a gluino has a distinct property in the collider experiments; a gluino 
with a decay length of CTg > 0{1) mm leaves a visible displaced vertex (DV) in the detectors, which 
greatly helps the gluino search. At present, however, there have been no dedicated searches from 
this aspect so far.^ 

The ATLAS collaboration has searched for DVs in the region of \z\ < 30 cm and r < 30 cm 

similar discussion has been recently given in Ref. [56] based on the CMS displaced dijets results [57], though 
their constraint is much weaker than ours. As we will discuss below, the ATLAS DV search [46] exploits the missing 
energy trigger, while the CMS search does not. In addition, the CMS dijet search requires large scalar sum of jet 
transverse momenta, which is not effective when the mass difference of bino and gluino is small. For these reasons, 
at present, the ATLAS search offers better sensitivities than the CMS one. 
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Figure 4: Current constraints (red and solid lines) and future prospects (blue and dashed lines) for 
the gluino searches. Favored region for the DM relic abundance is also shown in black lines for the 
cases of m = 50, 100, 200, and 300 TeV, with AM chosen so that the thermal relic abundance of the 
bino equals to the current observed DM density. We also show the current constraint [6] and future 
prospect of the 14 TeV LHC run [55] from the search for the prompt-decay gluino in horizontal red 
solid and blue dashed lines, respectively. 

in the inner detector [46, 58, 59], where z-axis points along the LHC beam line and r denotes the 
radial coordinate in the plane perpendicular to the ^-axis. They use the DVs reconstructed only in 
the air-gap region, namely, discard the DVs reconstructed within the material layers. This leads to 
signihcant background reduction. The signal region for the DVs is dehned such that the number of 
tracks associated with the DV is larger than four and ruDv > 10 GeV, where ruDv is the invariant 
mass of the tracks evaluated with the charged-pion mass hypothesis. Since they have observed no 
event in the signal region, they have given an upper limit on the long-lived gluino production cross 
section, which is interpreted as bound on the gluino mass in the high-scale SUSY scenario with a 
hxed neutralino mass of 100 GeV [46]. 

We re-interpret this low tudv search result in the case of the degenerate bino-gluino system, and 
obtain constraints on the bino-gluino coannihilation scenario, which is shown in Fig. 4. Here, the 
red and blue bands (from cTg = 1 mm to 1 m) show the estimated sensitivities of the DV search with 
the total luminosity of 20 fb~^ at the 8 TeV running and with 300 fb~^ at 14 TeV, respectively. The 
upper lines of these bands are for the cases where only the trigger efficiency is taken into account, 
which are simulated with HERWIG 6 [60] and AcerDET [61] to be 40% for 8 TeV with the 
threshold of the missing energy of 100 GeV, and 15% for 14 TeV with the missing energy trigger of 
200 GeV. Their dependence on the mass of gluino is only a few percent level. The lower lines, on the 
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other hand, correspond to the reconstrnction efficiency for DVs that is estimated from Refs. [58], 
where the long-lived nentralino decaying to two qnarks and one mnon is discussed in the i?-parity 
violating SUSY scenario. The reconstruction efficiency for the 108 GeV nentralino is about 20% 
of that for the 494 GeV nentralino in this case; we use this 20% for the lower lines, which gives 
conservative limits rather than the previous ones. In Fig. 4, we also show the favored region in terms 
of the DM relic abundance in black lines for the cases of fh = 50, 100, 200, and 300 TeV. Here, the 
bino-gluino mass difference AM is taken such that the thermal relic of bino DM explains the correct 
DM density. This reads that the present LHG data have already constrained a considerable range 
of parameter region consistent with the bino-gluino coannihilation scenario. This constraint is in 
fact much stronger than the ordinary limit from the searches of promptly decaying gluinos, which 
are based on only jets and missing energy [6, 7]. This constraint is indicated by the red and solid 
horizontal line in this hgure. The 14 TeV LHG running can further probe this scenario and reach 
Ms 2.5 TeV when cTg = 0(1 — 10) cm; this sensitivity is better than that by search with only jets 
and missing energy [55] (shown in the horizontal blue dashed line in the above hgure) by almost a 
factor of two. 

In addition, the ATLAS collaboration searches for massive charged meta-stable particles, such 
as i?-hadrons [62], with an another approach [47]. A characteristic feature of such particles is that 
they are produced with relatively low velocities, f3 = v/c < 1. This signature can be seen by means 
of large energy loss, dE/dx, in the ATLAS Pixel detector. Here, we note that this analysis requires 
gluinos to form charged i?-hadrons. Although the estimation of the charged hadronization fraction 
of gluinos may suffer from large theoretical uncertainty, this search offers the best sensitivity for 
CTg > 1 m. In Ref. [47], the result of this search is given as limits on the gluino mass in the case of 
AM = 100 GeV. We use the trigger efficiency given there for our computation for the 8 TeV case, 
and estimate the efficiency for the 14 TeV case by re-scaling it with a factor obtained by simulations. 
The red and blue solid curves in Fig. 4 show the estimated sensitivities of this search with 20 fb~^ at 
8 TeV and with 300 fb~^ at 14 TeV, respectively. We hnd that the searches of heavy stable charged 
particles give the most stringent constraints when cTg > 1 m, and are complementary to the DV 
searches. In particular, they are of importance when the scalar mass scale is relatively higher, say, 
a few hundred TeV. 

5 Conclusion and discussion 

In this paper, we study the bino-gluino coannihilation in the high-scale SUSY scenario. We have 
found that the squark mass scale cannot be too large for the coannihilation to work well. The upper 
bound on the squark mass is 200-1000 TeV for the gluino mass 1-8 TeV. Actually this mass scale is 
coincident with the prediction of the spectrum often called the spread or mini-split SUSY [13-21]. 
This constraint will provide a new perspective on the model-building to realize such mass spectrum. 

We also discuss the LHG signatures of this scenario. Because of the small mass difference between 
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the bino LSP and gluino, which is necessary for coannihilation, as well as heavy sqnark masses, the 
glnino decay length is considerably prolonged. Despite the small jets and missing energy activity, 
the DV and i?-hadron searches can efficiently probe snch long-lived glninos. If the sqnark mass scale 
is higher than abont 100 TeV, the cnrrent lower bonnd on the glnino mass is aronnd 1.2 TeV. The 
13/14 TeV LHC 300 fb“^ stage is expected to be able to explore glninos with a mass of ~ 2 TeV. 

Let us speculate possible sensitivities for much higher energy machines. For glninos with the 
decay length longer than 0{1) mm, a mass of 4.5 (10) TeV can be probed using a ^/s = 33 (100) TeV 
running proton collider with the integrated luminosity of 100 fb“^, provided that the background 
is sufficiently small and the detection efficiency of glninos is the same as that of the current LHC 
detector. This estimation may, of course, be too naive. Further detailed studies should be dedicated 
to see more precise prospects for such colliders, though we expect that they can probe the most of 
parameter space of the bino-gluino coannihilation. 

Lastly, we discuss the possibility of other gaugino coannihilation scenarios. As in the case of 
the current study, the small mass difference and heavier sfermion scale easily make the next LSP 
live long. For instance, in the case of the wino and gluino coannihilation, we may observe very 
exotic signatures; if the gluino lifetime is long enough, the gluino can carry the charged wino to the 
LHC trackers. In this case, we may observe displaced and disappearing tracks of the charged wino. 
The large gluino production cross section and the long-lived nature of the charged wino make it 
rather easy to look for this scenario in the LHC experiments. Another very interesting and plausible 
possibility is bino-wino coannihilation. This spectrum can be relatively easily realized even in the 
minimal anomaly mediation model. In this case, we may have another long-lived particle, which 
may play an important role at the LHC searches. A detailed analysis for this scenario will be done 
elsewhere [63]. 
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